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FOREWOED 


This  study  represents  the  joint  effort  of  Dr.  Gerhard;  W.  Braun, 
Chief  Scientist,  Deputy  for  Development  and  Test,  Air  Force  Missile 
Development  Center,  and  Dr.  Harald  A.  Melkus,  Senior  Aerospeuse  Engineer, 
Fluid  and  Flight  Mechanics,  of  the  Directorate  of  Research  Analyses, 

Air  Force  Office  of  Scientific  Research,  both  at  Holloman  Air  Force 
Beise,  Nev  Mexico.  The  analysis  vas  conducted  in  support  of  Project 
6876,  Track  Facility  Development,  assigned  to  the  Directorate  of 
Guidance  Test  Facilities. 

The  results  of  the  analysis  were  presented  at  the  UNESCO  Symposium 
on  Nonlinear  Oscillations,  held  in  Kiev,  USSR,  in  September  196I. 

Sensitive  guidance  con^nents  and  subsystems  are  affected  by 
vibrations  of  the  rocket  sleds  used  in  captive  flight  testing.  A 
quantitative  analysis  of  the  vibrations  induced  during  the  motion 
of  a  sled  along  the  rail  is  a  prerequisite  for  finding  mestas  of  reducing 
undesirable  vibrations  of  rocket  sleds.  Answers  must  be  obtained 
regarding  the  Influence  of  damping,  slipper  pl^  and  aerodynamic 
stabilization  on  the  dynamic  behavior  of  the  sled. 


ill 


ABSTRACT 


Some  resiilts  are  presented  of  an  analog  computer  study  concerned 
with  the  motion  of  a  rocket  sled  configuration  within  a  prescribed 
velocity  profile.  The  motion  of  the  sled  was  sidewise  restricted  by 
slippers  running  on  rails. 

The  dynamic  properties  of  the  sled  were  studied  under  Idealized 
conditions;  it  was  assiuned  that  the  moving  body  has  only  two  degrees 
of  freedom;  namely,  lateral  displacement  and  yaw.  The  body's  weather¬ 
cock  stability  was  varied  from  negative  to  positive  values  as  it  was 
assumed  that  this  aerodynamic  property  influences  the  mode  of  motion. 
The  effects  of  damping  and  that  of  play  between  slipper  and  rail  were 
studied  to  obtain  indications  as  to  how  these  parameters  affect  the 
motion. 


This  report  is  approved  for  publication. 


AMES  H.  RITTER 
blonel,  USAF 

Director,  Research  Analyses  Directorate 
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NCMENCLATURE 


coordimtes  of  the  c.g.  of  the  sled 
heading  angle  of  sled 
time,  seconds 

rail  gage  D  =  7  ft  (2.134  m) 

thrust,  if  positive,  or  brahing  force,  if  negative,  lb 

cross  forces  on  the  slippers,  positive  in  direction 
of  positive  y-axis,  lb 

magnitude  of  velocity  of  the  center  of  gravity  of  the 
sled  in  x-direction,  ft  sec-^ 

density  of  air,  lb  ft"'*  sec® 

aerodynamic  cross  force,  lb 

moment  of  inertia  of  the  sled  around  the  vertical  axis 
throxigh  the  center  of  gravity,  lb  ft  sec® 

spring  constant,  lb  ft"^ 

angle  of  yaw  of  the  sled  (p  =  ^-t^) 

displacement  of  slippermark  in  direction  of  positive 
y,  ft 

inner  and 


Subscript 

Subscript 

Subscript 

Subscript 

Subscript 

Subscript 


ou-ter  play  of  slippers,  ft 


Subscripts 


1  refers  to  left  front  slipper 

2  refers  to  right  front  slipper 

3  refers  to  left  rear  slipper 

4  refers  to  right  rear  slipper 
P  refers  to  front  beam 

R  refers  to  rear  beam 


ON  THE  NONUNEAE  YAW  MOTION  OP  A  ROCKET  SLED 


I.  INTRODUCTION 

In  missllry,  environmental  tests  have  become  of  ever-increasing 
iii5)ortance .  It  is  well  known  that  systematic  errors  are  generated  in 
high  precision  Instnunents— such  as  guidance  systems—when  subjected 
to  the  shock  and  vibration  environment  of  a  missile.  Their  con^iensai- 
tion  requires  a  calibration  under  closely  simulated  conditions  or,  in 
other  wcmds,  a  proper  reproduction  of  the  missile  environment. 

In  search  for  means  to  simulate  the  missile  environment  the 
similarity  between  a  supersonic  sled  and  a  missile  has  been  early 
recognized.  As  a  matter  of  fact  one  may  call  a  rocket  sled  a  rail 
guided  missile.  Figure  1  and  Figure  2  are  typical  exan5)les  of  rocket 
sleds  and  Illustrate  this  similarity. 

There  are,  however,  in^wrtant  differences  between  missiles  and 
sleds  which  cause  also  differences  in  their  internal  environments. 

The  most  conspicuous  differences  are  perhaps  the  presence  of  the 
slipper  beams  and  slippers  as  well  as  the  absence  of  tailplanes  for 
the  sled.  !Pherefore,  one  can  expect  that  the  main  difference  in  the 
environment  of  the  two  vehicles  is  introduced  by  forces  generated  in 
these  two  elements. 

For  a  quantitative  appraisal  of  effects  caused  by  these  forces, 
the  problem  was  sinplified  in  the  following  manner:  first,  it  was 
assumed  that  the  sled  moves  on  straight  rails,  and,  second,  its 
motion  is  limited  to  two  degrees  of  freedom--that  of  a  sideward  trans¬ 
lation  and  that  of  a  rotation  about  a  vertical  axis.  The  case  of 
straight  rails  is  of  practical  interest  as  it  represents  the  ideal 
condition  one  might  closely  approach  by  a  most  careful  aligning. 


grinding  and  smoothing  the  rails  as  has  been  done  for  quite  some  time 
on  the  Holloman  track.  At  the  same  time#  it  appears  that  the  confine¬ 
ment  to  straight  rails  does  not  represent  a  serious  handle apj  it  will 
be  seen  that  certain  typical  characteristics  of  the  motion  would  occur 
even  if  the  rails  were  perfectly  straight.  With  regard  to  the  other 
simplification  the  circumstances  are  similar.  Here  again  it  can  be 
expected  that  the  typical  properties  of  the  motlon—important  for  the 
detection  of  parameters  affecting  the  stability  of  a  sled— will  be 
preserved  to  a  great  extent. 

In  search  for  means  to  improve  the  stability  the  aerodynamic 
properties  of  the  moving  sled  were  incoiporated  in  the  analysis.  As 
a  rule  a  sled  without  a  tailplane  has  no  weathercock  stability;  in 
order  to  keep  a  sled  on  a  straight  path,  rather  high  and  frequent 
lateral  slipper  forces  are  required.  On  the  other  hand,  by  equipping 
a  sled  with  fins  one  might  expect  that  the  generated  aerodynamic  force 
will  contribute  to  the  stabilization  of  its  motion.  The  behavior  of 
a  sled  should  become  similar  to  that  of  an  aerodynamic  ally  stabilized 
missile  in  free  flight  provided:  (l)  the  sled  is  equipped  with 
sufficiently  large  fins,  and  (2)  the  pl^  between  slippers  and  rails 
is  made  so  large  that  lateral  disturbances  are  not  strong  enough  to 
cause  the  slippers  to  touch  the  rails  from  the  side. 

For  the  reeisons  stated  above  the  ideeaized  case  of  a  sled  with 
different  degrees  of  weathercock  stability  moving  on  straight  rails 
and  confined  to  small  translations  in  lateral  direction  as  well,  as  to 
small  yaw  angles  became  the  subject  of  an  investigation. 
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FIGURE  1.  DUAL  RAIL  GT.TiT> 


FIGURE  2.  MONORAIL  SLED 
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II.  DESCRIPTION  OF  THE  SLED  AND  RAILS 


For  the  environmental  tests  mentioned  before  a  sled  configuration 
consisting  of  a  pusher  and  forebody  is  mostly  used.  Figure  3  is  a 
photograph  of  such  a  sled  configuration.  The  test  specimen  is  mounted 
in  the  forebody  and  the  only  connection  between  forebody  and  pusher  is 
the  coupling  point.  The  coupling  has  the  form  of  a  universal  Joint 
which  also  provides  a  lateral  degree  of  freedom  (Fig.  4).  Due  to  the 
ever-present  vibrations,  friction  forces  in  lateral  direction  can  be 
expected  to  be  small.  In  the  following  investigation  they  are  con>. 
pletely  neglected.  This  is  to  say  that  the  thrust  force,  transmitted 
from  the  pusher  to  the  forebody  is  assumed  to  be  always  aligned  in  a 
direction  parallel  to  the  rails. 

On  a  sled,  vibrations  are  produced  by  thrust  irregularities  and 
by  aerodynamic  effects  such  as  shedding  of  vortices.  Sled  vibrations, 
however,  are  increased  by  the  disturbances  produced  by  the  rails  and, 
during  the  braking  period,  by  the  braking  device.  Therefore,  much 
effort  is  being  spent  to  strai^ten  and  smooth  the  rail  to  an  extreme 
degree  of  perfection.  The  master  rail,  i.e.,  the  western  rail,  is 
aligned  with  an  accuracy  of  .005  inch  (.127  mm)  tolerance;  the  eastern 
rail  with  .01  inch  (.254  mm).  But,  since  absolute  strai^tness  will 
never  be  obtained,  and  since  lypical  disturbance  patterns  become 
developed  even  on  straight  rails,  the  problem  still  remains  to  provide 
the  sled  with  such  dynamic  characteristics  that  it  has  a  minimum 
response  to  rail- or  brake-induced  disturbances.  In  the  following 
analysis  the  rails  are  considered  to  be  perfectly  straight. 
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FIGURE  5.  DIMENSIONS  AND  PHYSICAL  CONSTANTS  OF  THE  INVESTIGATED  SLED  CONFIGURATION 


III.  SLED  DYNAMICS 

,  .  .1 

This  study  is  primarily  concerned  with  the  dynamic  "behavior  of 
the  forehody  only.  In  the  following,  the  word  "sled”  refers  to  this 
portion  of  the  sled  team. 

The  sled  configuration  selected  for  the  analysis  consists  of  a 
fin-stabilized  axisymmetric  body  suppoirted  by  four  elastic  beams  (Fig.  5)* 
At  the  end  each  beam  is  equipped  with  a  slipper  used  to  guide  the  sled 
on  the  rails. 

The  forces  considered  instrumental  in  the  description  of  the  sled 
motion  are  as  follows: 

1.  The  thrust  surplus  T*  =  T  -  Dg,  exerted  by  the  pusher  in 
a  direction  parallel  to  the  rails. 

2.  The  rail  Induced  cross  forces  Si,  .  .  .  S4  acting  in  a 
direction  normal  to  the  rails. 

3.  The  aerodynamic  side  force  Y,  conq>osed  of  forces  acting 
on  the  eocisymmetrlc  portion  of  the  sled  and  on  the  vertical  tailplanes. 
Jb^licitly,  a  drag  force  was  not  Introduced.  It  was  assumed  that  this 
conqwnent  of  the  pressure  and  friction  for*ces  always  passes  through 
the  center  of  gravity  of  the  sled  and  is  baleuiced  at  every  instant  by 

a  corresponding  fraction  of  the  thrust  force.  The  thrust,  in  turn,  is 
eissumed  adjustable  in  such  a  way  as  to  neutralize  the  drag  at  each 


(1)  This  is  an  approximation.  Actually,  the  thrust  force  has  the 
direction  of  the  acceleration  of  the  coupling  point  between  ths  pusher 
and  the  forebody,  so  that  the  consideration  of  a  coupled  system  of 
pusher  and  forebody  would  be  more  appropriate  within  the  frame  of  a 
con5>rehensive  study.  For  this  simplified  treatment  of  the  motion, 
however,  the  asstm^tlon  that  the  thrust  vector  is  always  parallel  to 
the  rails  appeeurs  sufficient  and  more  revealing. 
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instant  and,  in  addition,  to  generate  a  specified  velocity  profile  with 
constant  slopes  of  +8,  0  and  -8  g's  (Fig.  6).  This  velocity  profile 
deviates  from  that  of  a  standard  sled  run  (Fig.  7)j  on-  the  other  hand, 
by  confining  the  motion  to  specified  kinematic  conditions,  the  analysis 
became  simplified.  All  forces  considered  are  confined  to  the  plane  of 
the  four  slippers  which  also  passes  throu^  the  center  of  gravity  of 
the  sled.  This  simplification  eliminates  roll  effects,  often  observable 
in  practical  cases,  but  still  leaves  enough  important  and  interesting 
nonlinear  phenomena.  Finally,  it  was  assiuasd  that  the  run  occurs  on  a 
calm  day  and  there  is  no  feedback  of  the  x-forces  between  forebody  and 
pusher  sled.  The  balance  of  forces  in  y-direction  yields; 


my  =  ^  Fy 

=  Y  +  Si  +  S2  +  S3  +  S4 


-Yp  2 


^  +  Si  +  So  +  So  +  s. 


(1) 


In  this  equation  the  following  notations  are  used: 


m 

y 

p 

u  =  u(t) 

D 

P 

°y 

yp 

Y 

Sa#  83^  S4 


mass  of  the  sled 

displacement  of  the  c.g.  of  sled  in  y-direction 
air  density 

velocity  of  the  sled  in  direction  of  positive 
x-axis  (along  the  rails) 

rail  gage 

angle  of  yaw  of  the  sled 

derivative  of  the  side  force  coefficient,  Cy.  of 
the  sled  with  respect  to  p 

aerodynamic  cross  force 

cross  forces  on  the  slippers,  positive  in  direction 
of  positive  y-axis 
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VELOCITY  PROFILE  SIMULATED  ON  THE  COMPUTER  FOR  THE  STUDY  OF 
THE  SLED  MOTION 
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FIGURE  7.  VELOCITY  PROFILE  FOR  A  OIVEH  PUSHER  -  FOREBODI  CONFIQURATION 
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It  is  assumed,  in  general,  that  the  aerodynamic  force,  Y,  is 
proportional  to  the  dynamic  pressure 


<1 


f 


and  is  a  linear  function  of  the  angle  of  yaw,  p  .  The  force  is  e3q)ressed 
in  terms  of  a  side  force  coefficient  Cy,  reflecting,  in  essence,  the 
effects  of  sled  geometry.  The  variation  of  Cy  with  Mach  and  Reynolds 
nuinber  is  disregarded.  The  equation 


c 


y 


(2) 


was  used  to  describe  the  linearized  side  force  effects  where  the  angle  of 
yaw,  p  ,  is  given  by  the  relation 


(3) 


In  equation  (3)  denotes  the  heading  angle  of  the  sled,  i.e.,  the 
angle  between  its  longitudinal  axis  and  the  x-orientation,  U  =  U(t)  is 
the  velocity  of  the  sled  in  x-direction,  and  y  is  the  cross-velocity 
of  the  sled. 

The  slipper  forces  Si  .  .  .  S4  depend  upon  the  elastic  properties 
of  the  slipper  beams  and  can  be  eapressed  in  terms  of  the  slope  K  of 
the  S  =  f(e)  curves,  displaying  the  relationship  between  the  side 
force  S  and  the  slipper  displacement  e  .  The  slopes  of  the  forces 
generated  on  the  front  slipper  beams  were  assumed  to  be  equal  for  the 
two  front  beams  but  different  from  those  of  the  rear  slipper  beams.  The 
slopes  (spring  constants)  were  designated  Kj.  and  K|^,  respectively. 
Figure  8  is  a  schematic  representation  of  these  forces. 
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rails. 
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THE  COORDINATE  SYSTEM  AND  THE  SYSTEM  OF  FORCES 
ACTING  ON  THE  SLED 
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$4#  when 


The  nonlinear  variation  of  the  slipper  forces  Si 
plotted  as  a  function  of  the  displacement,  €,  is  caused  by  the  play,  e, 
which  must  be  provided  to  insure  a  proper  gliding  of  the  sled  on  the  rails 
Within  the  frame  of  this  study  it  was  assumed  that  the  total  play,  > 
has  an  inner  part,  >  of  equal  magnitude  for  all  four  slippers  eind  an 
outer  part,  €q  (€q  *  €;j^),  also  the  same  on  all  four  slippers.  When  e^' 
and  the  e^'s  for  all  slippers  are  respectively  equal,  the  c.g.  of  the 
sled  is  exactly  on  the  middle  of  the  rails  and  the  sled  longitudinal  axis 
coincides  with  the  abscissa  of  the  x  -  y  coordinate  system. 

The  displacement  of  the  slippers,  e,  can  be  expressed  In  terms  of 
the  displacement,  y  ,of  the  c.g.  of  the  sled  and  the  heading  angle,  . 
Inspection  of  Figure  9  yields  the  relations; 

fii  =  y  +  fi  sin  i|f  -  ^  (l  -  cos  i|f) 

€2  =  y  +  ii  sin  ^  (1  -  cos 

€3  =  y  -  ^2  sin  i|f  -  I  (1  -  cos  \jr) 

€4  =  y  ~  ^2  sin  Tjr  +  ^  (1  -  cos  i|f)  (4) 

Considering  that  the  heading  angle  is  a  small  quantity 

^%iax  rad)  equation  (4)  can  be  linearized.  With  sin  ijr  s 

and  cos  ^  s  1,  one  obtains:  • 


€1  -  €2  3  y  +  f i  y 

€3  =  £4  s  y  -  ^2  y  (5) 
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FRONT  SLIPPER  BEAM 


t 

S  NOTE:  The  slopes  Kp  of  the 
S  =  f(e)  curves  are  equal  at 
corresponding  points. 


TOTAL 


PLAT 


The  yawing  motion  of  the  sled  is  described  by  the  moment  equation: 


i  =  N 

=  ^  +  Si(ii  cos  t  -  I  sin  i|r) 

+  SgCii  cos  t  if) 

-  83(^2  cos  T|f-f  ^  sin  i|r) 
“  84(^2  cos  t  “  §  t) 
+  T*  ^3  sin  ijf 


(6) 


The  first  term  on  the  right  side  denotes  the  aerodynamic  moment 
induced  by  the  lateral  and  angular  displacements  of  the  sled  as  well  as 
by  its  rotational  velocity;  I_  and  T*  denote  the  sled  moment  of 
inertia  about  the  z  axis  and  the  thrust  surplus,  respectively.  ^  is 
the  angular  acceleration  of  the  sled  and  ^3  is  the  distance  between 
the  coupling  point  and  the  c.g.  of  the  sled  (Fig.  3)» 

!Ehe  slipper  forces  are  nonlinear  functions  of  ei,  62,  €3  and  €4, 
respectively.  They  are  illustrated  by  Figure  9«  The  forces  are  multiplied 
by  the  nonlinear  functions  in  the  parentheses.  Since  the  yaw  angle  stays 
between  small  limits  it  was  considered  sufficiently  ewjcurate  to  replace 
cos  y  by  one  and  to  neglect  ^  sin  y  in  coinparison  to  and 

This  reduces  the  equation  (6)  to  the  siB5)ler  form; 

y  =  I  U^3  +  (Si  +  82)  ii  -  (S3  +  84)  ia  f  T*  i3  y  (6a) 
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For  small  angles  the  aerodynamic  moment  of  the  sled  can  be  thought 
of  as  the  sum  of  (l)  the  moment  generated  by  the  side  force  acting  upon 
the  axisymmetric  body  and  (2)  the  moment  produced  by  side  forces  on  the 
tail  planes  of  the  sled.  The  tail  force,  in  turn,  consists  of  two  parts; 
one  due  to  the  change  of  the  arigle  of  yaw,  the  other  due  to  the  rate  of 
rotation  of  the  sled.  We  combine  the  moment  generated  by  the  first  part 
of  the  tail  force  and  the  moment  induced  on  the  body  into  one  term 
(since  both  are  functions  of  the  aerodynamic  yaw  angle  p)  in  that  we 
write  for  this  portion  of  the  total  aerodynamic  moment 


Ni  =  cni  I 


Ab 

"d^ 


+ 


at  ^ 

^  Jt  D3 


P  •  I 


(7) 


In  this  equation  the  following  notations  are  used: 


the  slope  of  the  side  force  coefficient  generated 
by  the  body  of  the  sled. 


-  the  distance  between  the  center  of  pressure  and 
the  c.g.  of  the  sled,  positive  when  the  center  of 
pressure  is  located  aft  the  c.g. 


the  slope  of  the  side  force  coefficient  generated 
by  the  tailplanes 


-  distance  between  the  center  of  pressure  and  the 
c.g.  of  the  sled,  positive  when  the  center  of 
pressure  is  located  aft  the  c.g. 


reference  area  of  the  sled 
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area  of  the  tall  planes 


Aijx 

-  numerical  factor  describing  the  static  lateral 
^  stability  of  the  sled,  i.e.,  the  magnitude  of  the 

arestoring  moment  per  unit  angular  deflection  p 


The  rate  of  rotation  makes  a  similar  contribution  to  the  tail  force  as 
•an  apparent  increase  of  the  angle  of  yaw.  If  the  sled  rotates  with  the 
angular  velocity  y  and  the  distance  between  the  center  of  pressure  on 
the  tail  plane  and  the  center  of  gravity  of  the  sled  is  iq,,  the  lateral 
velocity  con5>onent  of  the  tail  is  equal  to  This  is  equivalent  to 

the  change  of  the  angle  of  yaw  of  the  amount  yjBij  U"^.  Thus,  the  rotation 
produces  a  force  on  the  tail  plane  idiich  is  proportional  to  the  angle 
yiT  ^  ,  and  we  can  write  the  corresponding  moment  about  the  center  of 
gravity  of  the  sled  in  the  form 


Na  =  Cnafu^DS 


Jlljt 

I>3 


yi^  I 


■Aqj 

D3 


y  I 


•  C  •  “*“ 

y  2 


y 


(8) 


where  c^  is  a  function  of  the  sled  velocity  U. 

Equation  (8)  describes  the  damping  effect  generated  by  the  tailplanes 
of  the  sled.  A  corresponding  aerodynamic  effect  of  the  anguleur  velocity, 
y  ,  induced  by  the  sled  body,  was  xieglected  because  of  its  smallness. 
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The  moment  equation  (6)  can  now  he  written  as  follows: 


^  - 


cnp  (^f  “  u  )  +  t 


I  U®d3 


+  (Si  +  Sz)  ii  -  (S3  +  S4)  ifi  +  T  is  i 


(6h) 


in  which  the  second  and  third  terms  on  the  right-hand  side  represent  the 
nonlinear  variation  of  the  moment  caused  hy  the  sled  displacement  y  and 
the  heading  angle  ijr  . 

Equation  (6b)  together  with  equation  (l)  constitutes  a  set  of 
simultaneous  •  differential  equations  of  second  order,  which  describes  the 
motion  of  the  sled  in  the  x-y  plane. 

IV.  APPROACH 

Even  for  the  idealized  cases  for  which  the  sled  motion  has  been 
confined,  an  analytical  treatment  of  the  problem  by  modem  analytical 
means,  known  to  the  authors,  appears  to  be  rather  involved.  Thus,  in 
order  to  obtain  a  quantitative  appraisal  of  the  various  parameters 
affecting  the  mode  of  motion  of  the  sled  and  its  stability,  several 
numerical  exan5)les  were  treated  by  simulating  the  sled  motion  on  a  PACE 
analog  con^juter.  For  practical  purposes  the  numerical  accuracy  of  an 
analog  conputer  with  an  estimated  error  of  not  more  than  .5  percent  was 
considered  by  far  sufficient. 

Furthermore,  the  stability  of  motion  was  examined  for  a  modified 
sled  configuration,  i.e.,  a  system  for  which  the  slipper  forces  can  be 
considered  linear  functions  of  displacement. 
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V.  ANALOG  COMPUTER  PROCEDURES 

Standard  techniques  were  used  to  mechanize  the  equations  of  motion 
on  the  analog  computer.  The  coiiiputer  time  was  chosen  to  he  eight  times 
real  time.  The  summing  aM  integrating  amplifiers  were  of  the  chopper- 
stabilized,  high  precision,  D.C.  type.  The  multipliers  were  all  electronic, 
using  the  time  division  principle. 

The  discontinuities  in  the  displacement-force  diagram  for  the 
front  and  the  rear  slippers  were  simulated  by  a  pair  of  smqjlifiers  with 
diodes  in  the  feedback  and  input  biaaes,  determining  the  cut-in  voltage 
of  the  outputs.  The  summing  amplifiers,  siamming  the  inputs  of  a  pair  of 
these  biased  diode  anplifiers,  will  give  an  output  voltage  only  if  the 
input  exceeds  the  bias  in  either  polarity.  Thus,  the  bias  voltage  deter¬ 
mines  the  ”dead”  space  which  corresponds  to  the  slipper  play  €•  Circuits 
of  this  type  eliminate  the  errors  due  to  contact  potential  of  the  diodes. 

.  A  difficulty  encountered  in  the  process  of  motion  simulation  con¬ 
sisted  in  finding  a  suitable  conpromlse  in  scaling  the  very  tight  spring- 
force  loops  and  the  relatively  loose  aerodynamic  loops,  ^e  time  scale 
of  eight  to  one  and  a  gain  of  approximately  400  for  the  spring-forces,  as 
they  ane  fed  into  the  integrators  for  y  and  p,  brought  a  satisfactoiy 
solution.  For  such  high  integrator  gains  a  certain  amo\nj.t  of  drift  is 
inevitable.  Physically,  this  can  be  interpreted  as  a  moiierate  side  wind. 

The  conputer  diagram  showing  the  simulation  is  given  in  Figure  10. 
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VI.  THE  LINEARIZED  SYSTEM 

For  the  linearized  system  the  slipper  forces  Si 


«  •  • 


S4  can  he 


following 

way  (see, 

e.g.. 

eq,uation  (5)  ): 

Si  =  - 

Kp 

€1  =  - 

Kp(y 

+ 

i) 

Sa  =  . 

Kp 

€2  =  - 

KpCy 

+  ii 

t) 

S3  =  - 

Kr 

€3  = 

KR(y 

-  £z 

t) 

84  =  - 

% 

€4  =  - 

KrCy 

“  £2 

t) 

(9) 


where  Kp  and  Kp  denote  the  modified  spring  constants. 

With 

^  5  ^  5  -  u  >  %(♦-  §  ) 

^n  P  •  •  ♦ 

N2  =  ”■? - =  N*  V  “  y  U® 

Sy  2  V  ? 

one  obtains  the  following  set  of  simultaneous  differential  equations 
of  second  order; 


m  y  +  Yp  Uy  +  2  K  y  +  (2  L  -  Yp  =  0 
\  t  -  U®  y  +  1^2  M  -  T*  ig  -  Np  U®  y  +  Np  U  y  +  2  L  y  =  0 


(10) 
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where 


K  =  Kp  +  Kj^ 

L  =  Kp  -  Kj^  iz 

M  =  Kp  +  Kj^  <^2^ 

the  coefficients  of  the  quartic 

A4  +  Ag  X®  +  Aa  >2  +  Ai  Tv  +  Ao  =  0  (ll) 

are  the  following  egression 
A4  =  ®  Iz 

Ag  =  U(Yp  Ig  -  m  Ny  U) 

Aa  =  m(2M  -  T*  ig)  +  2K  Ig  -  U^Cm  Np  +  Yp  U)  (12) 

Ai  =  U  j^Yp  (2M  -  1?*  ig)  -  2K  U  -  2K  Np 

Ao  =  2K(2M  -  T*  ig)  -  +  u2(2L  Yp  -  2K  Np)  . 

The  evaluation  of  the  coefficients  A4  .  .  .  Aq  for  various  values 

of  the  spring  constant,  K  =  Kp  +  Kp  ranging  from  K  =  6  X  10®  Ih  ft" 

to  about  K  .=  9  ^  10^  ft”^  has  shown  that  the  Routh  stability 
conditions 

A4J  Agj  Aa;  Ai;  Aq  >  0 

and  (13) 

Ai(Ag  Aa  —  A4  Ai)  —  Ag^  Aq  >  0 
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are  fulfilled.  Indicating  a  stable  system.  The  numerical  values  of  the 
various  A-coefficients,  the  resulting  natural  frequencies  of  the  motion 
and  the  damping  exponents  for  a  representative  configuration  vith 

Yp  =  0.051  lb  ft”^  sec^  (cyp  =  l.O) 
w  =  -0.02675  lb  ft“^  sec^  (cn-  =  -0.075) 

p  P 

=  -0.2283  It*  ft“^  sec”®  (Uc^  =  -0.64) 

are  con5>iled  in  the  table  on  page  27. 

VII.  RESULTS 

Some  of  the  results  obtained  in  the  form  of  coB5)uter  graphs  are 
displayed  in  Figures  11  through  I3,  It  should  be  noted  that,  for 
purposes  of  slii5)lificatlon,  the  assim^tion  of  equal  inner  and  outer 
play  coii5)onents  (ei  =  Cq  =  s)  was  iised  throughout  the  study.  Thus, 
whenever  the  front  or  the  rear  beam  of  the  sled  becomes  engaged  with 
the  rails  (because  of  a  lateral  displacement  or  a  rotation  of  the  sled) 
both  halves  of  the  pertinent  beam  were  equally  stressed.  It  is  possible, 
therefore,  to  substitute  for  the  spring  constants  of  the  two  halves 
of  the  front  beam,  a  •‘combined”  spring  constant,  Kj.  ,  and  for  those  of 
the  rear  beam  a  corresponding  constant,  Kj^. 

To  alleviate  coitqjarison,  six  con^uter  graphs,  each  of  them 
representing  a  complete  sled  run  of  about  27*4  seconds  duration,  are 
assembled  in  one  figure.  Each  graph  displays  eight  variables  which 
characterize  the  lateral  motion  of  the  sled  and  their  variation  with 
time,  t.  Plotted  are  from  top  to  bottom; 

1*  Yp  =  Yp  (t)  the  lateral  displacement  of  the  end  point  of 

the  front  beam,  if  rigid 

2.  yj^  =  yj^  (t)  the  lateral  displacement  of  the  end  point  of 

the  rear  beam,  if  rigid 
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3.  y  =s  y  (t)  the  lateral  displewement  of  the  c.g, 

of  the  sled 

4.  p  =  p(t)  the  yaw  angle  of  the  sled 

5*  y  =  y(t)  the  lateral  velocity  component  of  the 

c.g.  of  the  sled 

6.  •|r  =  ijf(t)  the  angular  velocity  of  the  sled 

7»  Sj,  =  Sp  (t)  the  force  acting  on  the  rear  end  of 

the  sled 

8.  Sp  =  Sp(t)  the  force  acting  on  the  front  end  of  the  sled 

The  time  scale  is  displayed  at  the  bottom  of  the  graphs,  each  unit 
representing  one-eighth  of  a  second.  (As  mentioned  before,  the  con5)uter 
time  was  chosen  eight  times  real  time;  thus  a  sled  run  of  2J,h  seconds 
duration  is  composed  of  27.4  X  8  =  219.2  units,  i.e.,  requires  about  220 
seconds  on  the  con^juter. ) 

For  a  period  of  about  6.22  seconds  (~  50  units)  the  sled  was  sub¬ 
jected  to  a  constant  acceleration  of  8  g's  which  produces  a  velocity 
increment  of  A  U  «  l600  ft  sec”^.  In  the  study,  all  sled  runs  started 
with  an  initial  velocity  of  Uq  =  200  ft  sec”^  (to  avoid  conqjutational 
difficulties  as  the  computer  had  to  evalxmte  expressions  containing 
the  sled  velocity  in  the  denominator)  so  that,  at  the  end  of  the  accelera¬ 
tion  period,  the  velocity  of  the  sled  had  reached  the  value  U  =  1800  ft  sec"^. 
(Fig.  6)  For  about  15  seconds  (120  units)  this  velocity  was  held  constant; 
then  by  applying  a  negative  acceleration  of  -8  g's  the  sled  was 
decelerated  until  its  velocity  had  reached  the  200  ft  sec”^  mark. 


The  upper  row  of  computer  graphs  in  Figure  11  displays  three  sled 
runs  in  which  the  slipper  play,  e,  was  varied  from  e  =  .0013  ft 
to  e  =  .0039  ft  while  the  aerodynamic  coefficients,  Cy  ,  c^  and  Ucl 
as  well  as  the  spring  constants  and  were  held 

constant  (2Kj.  =  Kjj  =  4  X  10®). 
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table: 


-1 

U  ft  sec 

T*  lb 

500 

9973 

1000 

13011 

1500 
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9114 
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o  in 

o 

1  r— 1 

"  4.  X 

&.• 

O  II 

II 

o  ->1^ 

•  o 

H  11 

oa  CM 

o 

A4  X  10"“*  lb2  sec^ 

A3  X  10  ■*  ib^  sec® 

Az  X  10“S  Ib^  sec® 

Ai  X  10""®  lb®  sec^ 

Aq  X  10"^®  lb® 

cUi  sec 

1  -1 

cjc^  sec 

cps 

rig  c:ps_ 

sec_ 

Vg  sec" 

2.958 

2.5648 
1.8655 
7.3012 
2.5428 
207.7 
141.2 
33.06 
22.47 
-.  2740 
-.1595 

2.958 

5.1295 

1.8655 

14.5928 

2.5356 

208.0 

l40.8 

33.10 

22.41 

-.5477 

-.3193 

2.958 

7.69^3 

1.8656 

21.8654- 

2.5236 

208.4 

140.2 

33.16 

22.31 

-.8208 

-.4798 

2.958 

7.69^3 

1.8676 

21.9083 

2.5304 

208.4 

i40.3 

33.17 

22.33 

-.8197 

-.4809 

LA 

O  ^ 

.*  s 

MO 

01  •  ^ 

^  ‘ 

u  II 

II 

o 

H  g  II 

" 

^  CM 

o 

A4  X  10"^  lb®  sec-* 
A3  X  10  ■*  lb®  sec® 
As  X  10"®  lb®  sec® 

Ai  X  10"“^  lb®  sec^ 

Ao  X  10"^^  lb® 

ciii  sec  ^ 

o:^  sec"^ 

cps 

Ha  cps^ 

Vi  sec_ 

Vg  sec"^ 

2.958 

2.5648 

1.8486 

7.1652 

2. 4616 

65.77 
43.86 
10.47 
6.98 
-. 2765 
-.1570 

2.958 

5.1295 

1.8489 

14.2352 

2.3897 

66.51 

42.73 

10.59 

6.80 

“.5503 

-.3167 

2.958 

7.6943 

1.8500 

21.1318 

2. 2697 
67.64 
40.95 

10.77 

6.52 

-.8207 

-.4798 

2.958 

7.6943 

1.8484 
21.5430 
2.3370 
66.97 
41.96 
10.66  . 
6.68 
-.8048 
-.4958 

lA 

C^ 

O  0) 

1  s 

■  4  X 

&■•  t 

II 

•  •V  •■  PS 

H  5  II 

"  ^ 

S  CM 

O 

A4  X  10"*  lb®  sec-* 

A3  X  10  ■*  lb®  sec® 

Az  X  lO""^  lb®  sec® 

Ai  X  10"®  lb®  sec^ 

Ao  X  10"®  lb® 

OJi  sec_^ 

flc^  sec  ^ 

rij^  cps 

xig  c:ps_ 

V,  sec 

1  -1 

V-  sec 

2 

2.958 

2.565 

1.667 

5.755 

1.649 

20.86 

11.32 

3.32 

1.80 

-.2976 

-.1358 

2.958 

5.130 

1.683 

10.556 

0.929 

22.49 

7.88 

3.58 

1.25 

-.5865 

-.2806 

2.958 

7.694 

1.692 

13.448 

-0.27 

Unstable 

2.958 

7.694 

1.887 

17.743 

0.4o4 

24.77 

4.69 

3.94 

.75 

-.8429 

-.4577 

FIGURE  11.  VARIATION  OF  THE  SLED  MOTION  PARAMETERS  WITH  TIME 
Spring  Constants;  2  Kp  =  K|^  =  4  X  10®  lb  ft”^ 
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FIGURE  13.  VARIATION  OF  OIIE  Sm)  MOTION  PARAMETERS  WITH  TIME 
Spring  Constants;  2  Kp,  =  Kj^  =  4  X  10®  lb  ft“^j  Cy^  =  .5;  Cjj^  =  -.075 
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It  can  be  noted  that  a  periodic  motion  becomes  established  as  soon 
as  the  sled  reaches  the  U  =  const  portion  of  the  run.  Nearly  all 
motion  variables  exhibit  high  upper  limits.  Comparing  the  graphs  from 
left  to  right,  it  can  be  seen  that  with  increasing  slipper  play  a  gradual 
increase  of  the  amplitude  of  the  displacement,  y^,  ,  takes  place. 
Rememibering  that  yp  represents  the  totaJL  displacement  of  the  front 
beam  (i.e.,  the  play,  e,  and  the  deflection  of  the  spring,  Sp  Kp^  ) 
its  magnitude  is  not  indicative  of  the  elastic  force,  Sp,  which  is 
proportional  to  the  difference  A  =  yp  -  e.  Thus,  the  recorded  gradual 
reduction  of  Sp  with  increasing  e  does  not  contradict  the  above 
described  behavior  of  the  displacement  yp. 

On  the  other  hand,  the  variation  of  the  forces,  Sr  f  exerted  on  the 
rear  beam,  with  e  is  dlfferentj  these  forces  first  increase  with  increas¬ 
ing  €,  but  they  become  smaller  as  e  is  further  increased.  The  examina¬ 
tion  of  the  graphs  shows  that  this  behavior  is  related  to  the  mode  of 
motion  of  the  sled.  In  the  first  two  graphs  the  motion  can,  perhaps,  be 
described  as  a  rotation  about  a  "point”  located  several  feet  behind  the 
rear  beam  of  the  sled,  the  distance  increasing  with  increasing  e.  In 
the  third  graph,  however,  the  "point"  of  rotation  is  shifted  forward  and 
is  now  lying  between  the  c.g.  and  the  rear  beam.  Thus,  the  instantaneous 

forces  exerted  on  the  front  and  the  rear  beams  which  in  the  first  two 

\ 

cases  had  equal  signs  are  now  sictlng  in  opposite  direction. 

Finally,  it  can  be  noted  that  with  increasing  slipper  play,  e, 
a  considerable  reduction  of  the  dominant  frequency  is  taking  place. 

(n  s  l4  cps  for  €  =  .0013  ft  as  compared  to  n  s  10. 5  for  e  =  .0039  ft) 

The  lower  row  of  Figure  11  displays  the  behavior  of  the  motion 
variables  for  a  constant  value  of  the  slipper  play  e  (e  =  .0026  ft). 

In  this  arrangement  the  coefficient  of  the  aerodynamic  moment,  Cjj^,  was 
varied  from  -.5  to  -.15  while  Cy^  and  Uc^  were  held  constant.  The 
meaning  of  this  is  that  the  center  of  pressure  of  the  aerodynamic  load 
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-,ras  shifted  from  a  position  at  a  distance  of  7  feet  aft  the  c.g. 
(hypothetical  case)  to  a  position  2.1  feet  aft  the  c.g.  of  the  sled. 

It  can  be  observed  that  with  decreasing  distance  a  gradual  reduction 
of  forces  and  corresponding  frequencies  t alee s place. 

Figure  12  displays  the  effects  of  damping  for  various  values  of 
the  slipper  play  e.  The  spring  constants  are  chosen  to  be  five  times 
larger  than  those  in  the  previous  case.  The  aerodynamic  coefficients 
selected  for  this  configuration  are  those  of  the  representative  sled 
shown  in  Figure  5  (cy^  =  .5;  c^^  =  -.075).  The  upper  row  of  graphs 
displays  the  motion  of  a  sled  for  which  the  dan5>ing  coefficient  is 
assumed  to  be  zero,  while  the  lower  row  exhibits  a  configuration  having 
a  moderate  daii5)ing  coefficient  (Uc^  =  -.  40  to  -.50).  The  benefits 
obtained  by  dan^jing  are  obvious.  It  is  interesting  to  note  the  distinct 
variation  of  forces  exerted  on  the  front  and  the  rear  beam  with  damping) 
while  for  Uc^  =  0,  the  forces  on  the  front  and  the  rear  beam  are 
nearly  equal,  a  moderate  danqplng  coefficient  of  about  -.4  to  -.5  causes 
the  spring  forces  on  the  rear  beam  practically  to  disappear. 

The  con5)arison  of  the  frequencies  of  the  sled  runs  in  the  upper 
row  with  those  on  the  lower  row,  shows  that  dancing  also  causes  a  sub¬ 
stantial  reduction  of  the  dominant  frequency,  n;  the  n-values  decrease 
from  about  7*8  cps  (upper  row)  to  about  3*5  (lower  row). 

Figure  13  displays  the  effects  of  danq)ing  for  a  sled  equipped  with 
softer  springs  (l^  =  2  X  10  lb  f t  =  4  X  10®) .  Here  again  the  benefits 

resiilting  from  a  properly  selected  dan5)ing  coefficient  are  recognizable. 

It  should  be  noted,  however,  that  a  sled  equipped  with  softer  slipper 
beams  is  more  susceptible  to  the  magnitude  of  the  slipper  play  e#  than 
a  sled  with  stiffen  beams.  While  for  €  =  .0013  Toot  the  slipper  forces 
on  the  front  beam  are  exceptionally  small  (Sp  =  135  l^i  their 

magnitude  increases  continually  with  increasing  slipper  play.  (Sp  =  315 
Sp  =  100  lb  for  e  =  .0039  Tt) 
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To  obtain  a  better  insight  into  the  details  of  the  sled  motion  the 
time  scale  in  Figure  l4  was  expanded.  The  two  conputer  graphs  shown 
represent  the  first  portion  of  the  sled  run,  i.e.,  the  acceleration 
phase  and  the  beginning  of  the  U  =  const  phase.  It  can  be  noted  that 
the  displacements  y^,  and  have  opposite  signs;  every  time  the 
front  beam  moves  upward  the  rear  beam  moves  downward.  The  displacement 
of  the  c.g.  of  the  sled  is  relatively  small.  The  motion  can  be 
described  as  a  rocking  of  the  sled  about  a  '‘point"  located  in  the 
neighborhood  of  the  c.g.  of  the  sled.  Actually,  the  "point”  is  a 
relatively  small  area  to  which  the  instantaneous  center  of  rotation  is 
confined  and  in  which  the  center  is  following  a  complicated  path.  However, 
in  the  first  approximation  the  area  can  be  regarded  as  a  "point"  about 
which  the  sled  is  oscillating. 

The  path  described  by  the  end  point  of  the  rear  beam  is  nearly,  a 
sinusoid  with  an  aiiplitude  which  is  smaller  than  the  slipper  play,  e. 

No  elastic  forces  are  exerted  on  the  rear  beam.  Thus,  the  sled  can  be 
regarded  as  a  rigid  body  with  elastic  forces  acting  on  its  front  end  in 
intervals  whenever  the  slipper  stops  are  touching  the  rails,  and  with  a 
periodic  aerodynamic  force  with  its  center  located  1.05  feet  aft  the  c.g. 

A  detailed  analysis  of  the  recorded  modes  of  motion  of  the  sled 
could  not  be  ceirried  out  as  time  was  pressing  and  only  a  very  limited 
number  of  conputer  graphs,  suitable  for  an  analysis,  was  available  to 
the  authors  at  the  time  of  this  writing. 

The  variation  of  the  slipper  forces  S  and  of  the  dominant 
frequency  n  with  slipper  play  e  is  shown  on  Figures  15  and  I6, 
respectively.  The  spring  constants  and  the  aerodynamic  coefficients 
were  held  constant  (see  Fig.  15  and  I6).  Details  of  the  sled  motion 
which  correspond  to  these  charts  are  presented  in  the. upper  row  of 
Figure  11.  It  can  be  noted  that  for  smaller  e's  a  gradueil  reduction 
of  the  slipper  forces  on  the  front  beam  and  an  increase  of  forces  on 
the  rear  beam  takes  pleice  with  increasing  e .  This  behavior  is 
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cheiracteristic  of  that  mode  of  sled  motion  for  which  the  center  of 
rotation  is  located  aft  the  rear  heam.  As  soon  as  the  motion  pattern 
changes  (the  center  of  rotation  heing  located  between  the  front  and 
the  rear  slipper  beam)  a  considerable  reduction  of  the  rear  slipper 
forces  can  be  observed. 

Extrapolating  the  measured  frequencies  (Fig.  l6)  toward  e  =  0 
one  finds  that  with  no  slipper  play  the  sled  will  oscillate  with  the 
,  frequency  n  of  about  15.9  cps.  A  frequency  of  15.91  cps  was  obtained 
analsrbically  by  assuming  that  a  rigid  body,  having  a  mass  and  an  inertia 
moment  corresponding  to  that  of  the  sled,  is  supported  by  two  springs 
with  constants  of  the  same  magnitude  as  those  of  the  sled.  The  observed 
shape  of  the  sled  motion  which  corresponds  to  the  frequency  of  15.9  cps 
is  the  second  natural  mode  of  the  oscillating  body  (lower  frequency, 
center  of  rotation  located  aft  the  c.g. ). 

Figures  17^18  and.  19  show  the  effects  of  damping  for  various  values 
of  the  slipper  play  e  and  of  the  spring  constant  K.  It  can  be  seen 
that  in  practically  all  cases  a  considerable  reduction  of  the  slipper 
forces  on  the  front  and  on  the  rear  beam  taJses  place.  It  should  be  kept 
in  mind,  however,  that  the  true  aerodynamic  coefficients  may  be  sub¬ 
stantially  lower  than  the  values  assumed  in  this  study.  In  the  p-near-zero 
range  the  bomdary  layer  might  have  an  adverse  effect  on  the  sensitiveness 
of  the  sled  to  small  angular  deflections,  thus  reducing  the  magnitude  of 
the  aerodynamic  coefficients. 

Figure  20  shows  the  variation  of  the  slipper  forces  S  with 
spring  stiffness  K  for  various  values  of  the  slipper  play  e.  The 
mode  of  motion  of  the  sled  was  that  for  which  the  center  of  rotation  is 
located  between  the  slipper  beams.  It  appears  that  by  proper  adjustment 
of  the  slipper  play  e  and  of  the  spring  stiffness  K,  the  slipper 
forces  can  be  minimized. 
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Figure  21  shows  the  variation  of  the  dominant  frequency,  n,  with 
the  spring  stiffness,  K,  the  slipper  play,  € ,  helng  the  parameter.  It 
can  be  noted  that  larger  frequency  changes  take  place  only  in  the  range 
of  lower  spring  stiffness,  K.  At  higher  K's  only  minor  changes  in 
n  can  be  observed.  The  mode  of  motion  of  the  sled  is  the  same  as  that 
of  Figure  20.  For  comparison  purposes  the  variation  of  the  frequencies 
of  the  undamped  sled  configuration  was  also  included. 

VIII.  CONCLUSIONS 

The  results  presented  indicate  that  the  simulation  of  a  rocket 
sled  motion  on  an  analog  computer  exhibits  the  possibility  of  being 
efficient  in  the  attempt  to  understand  the  mechanism  of  the  sled  motion. 
Through  the  use  of  this  technique  some  tedious  computational  difficulties 
could  be  avoided  and  the  analysis  of  the  dynamic  properties  of  a  sled 
simplified. 

The  method  provides  answers  on  the  influence  of  various  parameters 
on  the  behavior  of  the  sled,  thus  permitting  conclusions  regarding  their 
optimization. 

The  data  presented  herein  must  be  considered  preliminary  as  only  a 
limited  number  of  cases  were  examined. 
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FIGURiO  IS.  VARIATION  OF  THE  SLIPPER  FORCES  S  '/ITH  SLIPPER  PLAY  e 

PARATtSTER:  SPRING  CONST.AWT  K 
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FIGURE  19.  VARIATION  OF  THE  SLIPPER  FORCES  S  '-aTH  SLIPPER  FLAY  e 
PARAMETER:  DAMPING  COPJFFICIENT  U  ; 
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